ERROR BOUNDS FOR QUASI-MONTE CARLO INTEGRATION FOR 
^eo WITH UNIFORM POINT SETS 

— . , SU HU AND YAN LI 

(N 

^ ' Abstract. Niederreiter Ul established new bounds for quasi-Monte Carlo in- 

O i tegration for nodes sets with a special kind of uniformity property. Let (X, £/, fi) 

^^ ' be an arbitrary probability space, i.e., X is an arbitrary nonempty set, si/ a cr- 

algebra of subsets of X, and // a probability measure defined on £/. The func- 
tions considered in HI are bounded /i-integrable functions on X. In this note, 
we extend some of his results for bounded //-integrable functions to essentially 
bounded ^-measurable functions. So Niederreiter's bounds can be used in a 
more general setting. 



o 






1. Introduction 



Let (X, £/, /J.) be an arbitrary probability space, i.e., X is an arbitrary nonempty 
>• ' set, £/ a cr-algebra of subsets of X, and yu a probability measure defined on £/. 

^ I Niederreiter [IJ established new bounds for quasi-Monte Carlo integration for 

in ' 

in, ^ 

y^ , ered in [[T]| are bounded /i-integrable functions on X. In this note, we extend 

O I some of his results for bounded yu-integrable functions to essentially bounded 

.^/-measurable functions. 

"^ ■ 2. Main results 

c^ ■ Let =Sfoo(X, £/,//) be the set of all essentially bounded ^/-measurable func- 

tions on X, two functions being identified if they differ only on a //-null set. 
For any ^-measurable function g on X, ||^||oo(esssup^g;;-|^|) denotes the essential 
supremum of \g\ (see P.346 of 0). For an extended real-valued function /, we 
define /+ = max{/,0} and /" = -min{/,0}. Notice that /+ > 0,/" > 0, and 
/ = /^ - /" (see P. 164 of ^). For a given nonempty subset ^ of £/, let L^^ 
be linear subspace of ^oo(X,£/,j^) spanned by the constant function 1 and all 
characteristic functions ;\fM5 M e M . For any / e =Sfco(^, ^,/^), let D(/, Lj^) be 
the =Sfoo distance from / to ^^, that is , 

D(/,L^) = inf,ei.^||/-Z|U. 
The following definition can be found in P.285 of |[I1. 
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Definition 2.1. Let {X, s^,ijl) be an arbitrary probability space, let ^ be a nonempty 
subset of s^ . A point set l^ofN elements oiX is called (^, //)-uniform if 



N 



J^ATmC^J = A{M; ^) = iu(M)N, for all Me^. 



;=1 



Let {X X ... xX,f/x ... X ^,/i X ... X/j.) be the product measurable space 



N 



N 



(see P.379 of ^). We can view a point set ^ = {Xi,...,Xn} as a point in 
X X . . . X X and ^ Z^=i /(^n) as a //-variable function on X x . x X . Since / e 

^ooiX, £/,iu), we have ^ Zti fi^n) e ^o,( X x . . . xX , ,g/ x . . . x ^ ,jU x ... x /u). 

N N !^ 

Let 

"rf = {(Xi,...,Xa,) e X x . xX\^ = {Xi,...,X/v}isan(^,/z)-uniform point set}. 

N 

Since f{Xi,X2...,Xi^) = Yjf=i XiwiXn) is a measurable function on Xx...xX, from 
Definition 12. II and Lemma 1 1.9 of [2], if ^ is a countable nonempty subset of 
£/, then ^ is a measurable set. 

Theorem 2.2. Lef (X, s^, jj.) be an arbitrary probability space. Let ^ be a count- 
able nonempty subset of s^. Let 

^ = {{Xi, ...,Xi^) £ X X . . . xX \S^ = [Xi, ...,Xfj}is an(^,p)-uniform point set}. 

N 

Then for any f e <ifoo(X, £^,fJ.), we have 

1 ^ r 

esssup^Xu...,x„wATTYjf^^"'>~ I fdfi\<2D(f,L^). 

Proof. We extend Niederreiter's proof for Theorem 1 of [1 J to our case. For any 
M £ .M and any (^,/i)-uniform point set 3^ = \X\, ...,Xt^]. we have 

1 /. A(M\^) r 

^ 2^XM{Xn) = = l^(M) = I XMdp 

by the definition of an (^,//)-uniform point set. 

For any / G =Sf^ and any (^,ju) -uniform point set J^ = {Xi, ..., X^^}, we have 



r 1 '^ 

idfi=-yiiXn). 
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Thus for any / e ^co(X, £/,ij.) and any (^,yu)-uniform point set ^ = {Xi, ..., X^} 
we have 

r 

fdn 



- y (/ - i)(x„) + - y i(xj -if- Ddiu - idfi 



n=\ 

for all l€L^. 
So 



1 " r 



1 r 

1 ^ r 

<eSSSUP(^^ X„)e^<flT7y,(/-0(^n)l+ 1/ - ^1^/^ 

1 ^ r 

<esssup(;f^_x„)g.^lT7y(/-0(^n)l+ l/-/|^// 

^^ t^ J{xeX||/-;|>||/-/|U) 

r i/-/irf// 

J(A:eXllf-/l<llf-/IUI 



+ 

*{^eX||/-/|<||/-/|| 

for all l&L^. 
Since 

N 



{(Xu ...,Xn) e "^ I I^ Yjif - l)iXn)\ > 11/ - l\U 

n=l 

1 ^ 

c {(Xi, ...,X^) e X>i^^^2<J I I- Yj^f - l){Xn)\ > 11/ - l\U 



c u;^^i(X X ... X {X„ e X I |(/ - l){Xn)\ > 11/ - /|U}... X X\ 
we have 



1 
/0^^^({(Xi, ...,X^) e ^ I I- J](/ - l){Xn)\ > 11/ - l\U) 

1 " 

< /o^^^({(Xi, ...,X^) e X2i^^^x^ I I- ^(/ - /)(X„)| > 11/ - /|U}) 

/V 

< J]m2^^^(X X ... X {X„ G X I |(/ - l)(Xn)\ > 11/ - /lloo}... X X) 

N 

= J]K{Xn eXWif- /)(X„)| > 11/ - l\U) = 0, 



n=l 
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the last equality follows from Fubini's theorem (see P.384 of Q). 
From the definition of =^oo-norm, we have 

1 ^ 






Also from the definition of ^co-norm, we have 

MUg^II/-/|>II/-/||oo}) = 0, 



thus 



So 



r 1/ - Mm = 0. 

J{:c€A:||/-/|>||/-/|U) 

1 ^ r 

esssup(;^, ;f^)<..^|-y /(XJ- fd^\ 
<II/-/|L+ r \f-l\dM 

J{xeX\\f-l\<\\f-l\\„\ 

<2||/-/|U 



for all / e L, 



Let ^ = {Ml, ..., Ml,} be a finite nonempty subset of £/ such that Mi, ..., M^ 
are disjoint and U^^^jA/,- = X. If / e ^co{X, £/,n), then / 6 ^^{X, £/\m,,IJ-\m,) for 
any I < i < k. Let 

(^ .(^) = |-inf.eM/-, if M{^ 6 M,ir « > 0}) = 0; 
^ lesssup^g^/"^, otherwise, 

|inf.eM/\ if /^({-^ e Mj\f-(x) > 0}) = 0; 

^ 1-esssup^g^ /", otherwise, 

for 1 < j < k. Define 

S.Af) = maxi<y<^(G//) - gj(f)). 

Corollary 2.3. Let (X, £/, fi) be an arbitrary probability space. Let ^ = {Mi, ...,M^ 
be a finite nonempty subset of s^ such that Mi, ..., Mj, are disjoint and U*_,M; = 
X.Let 

'€ = {{Xi, ...,Xn) e X X ■ X X \^ = {Xi, ...,Xi^]is an{^ , ij)-uniform point set). 

N 

Then for any f e ^co{X, £^ ,p.), we have 

1 ^ r 

esssup(x^_x^)e<AT-ry]f^^"^~ I fdlA^S^{f). 
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Proof. We extend Niederreiter's proof for Corollary 1 of [JJ to our case. Let 

for \< i<k, let 

k 

Since 

I, e M, I |/(,)-C,| > M>ZM)| c I, . M,W) > G,(/)|UU e M, | /(,) < ,//),, 

we have 

Ml. . M, I |/W-C,| > 5^<^^|) . Ml' e M, I /(,) > G,(/))).MI' . M, I /(,) < «,(/),) = 

by the definition of =Sfoo-norm. 
For t e Mj, we have 

esssup,g^.|/(0-/(OI 

= esssup,g^J/(0 - Cj\ 

<\(Gj{f)-gj{f)) 



Therefore 



Thus 



\\f-l\L<y.Af). 



Dif,L^)<^S,Af)- 



From Theorem l2.2[ we get our result. n 

Corollary 2.4. Let (X, ^, /j) be an arbitrary probability space, let ^ = {My , ..., Mj,] 
be a finite nonempty subset of s^ such that Mi,..., Mt are disjoint and u'^.^Mj = 
X.Let 

'^ = {{Xi,...,Xn) e Xx . xX \^ = {Xi, ...,XN]is an(^,p)-uniform point set]. 

N 

Then for any f e ^oo{X, s^ , /x) , we have 

^ N k 

^***"/'(x,,...,x„)6^#It7 y!/(^n) - fdp\ < VpiMjXGjif) - gj{f)). 
^ t^ Jx ^ 

Proof. We extend Niederreiter's proof for Theorem 2 of |[I1 to our case. From 
the definition of ^oo-norm, we have 

KMj)gjif)< f fdn<ix{Mj)Gj{f). 
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From the definition of uniform point set, we have 



1 ^ 



n=\ 



c u;)Li(X X ... X {Xn e Mj I gjif) > /(X„)}... X X), 
for I < j <k. 

Thus 

1 ^ 
IJ X ... xM{(^i, -,^iv) e "^ I KMj)gj(f) >t;J] -^^^"^^^ 



w 



< J]// X ... xM^ X ... X {Xn e My I gjif) > /(X„)}... XX) 



n=l N 



= J]iu({Xn e Mj I gjif) > f(Xn)}) = 

for 1 < J < fc, by the definition of gj{f) and Fubini's theorem. Similarly, 

1 ^ 
// X ... xMK^i, .■•,^;v) e ^ I KMj)Gj(f) < T^ Xi ^^^"^^^ = ^- 






Thus 

1 

esss 






X„eMj 

for 1 < 7 < ^. 
Finally, from 



N 

fdix 






we get our result. 



Acknowledgement: This work was partially supported by Postdoctoral Sci- 
ence Foundation of China. The authors are grateful to the anonymous referee for 
his/her helpful comments. 

References 

[1] H.Niederreiter, Error bounds for quasi-Monte Carlo integration with uniform point 
sets, Journal of computational and applied mathematics 150 (2003), 283-292. 

[2] E.Hewitt and K.Stromberg, Real and abstract analysis, GTM25, Springer- Verlag, 
Berhn Heidelberg, 1965. 



ERROR BOUNDS FOR QUASI-MONTE CARLO INTEGRATION FOR .S*„ 

Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 
E-mail address: hus04@raails . tsinghua . edu . en 

Department of Mathematical Sciences, Tsinghua University, Beijing 100084, China 
E-mail address: liyan_00®mails . tsinghua . edu . en 



